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AN ANALYSIS OF ABLATION-SHIELD REQUIREMENTS FOR MANNED REENTRY 

VEHICLES 


"By Leonard Roberts 


SUMMARY 

The problem of sublimation of material and ac- 
cumulation of heat in an ablation shield is analyzed 
and the results are applied to the reentry of manned 
vehicles into the earth’s atmosphere. The param- 
eters which control the amount of sublimation and 
the temperature distribution within the ablation 
shield are determined- and presented in a manner 
useful for enyineeriny calculation . ft is shown that 
the total mass loss from the shield durintj reentry 
and the insulation requirements may he yiren very 
sim ply in terms of the maximum deceleration of the 
vehicle or the total reentry time. 

INTRODUCTION 

The successful return of a vehicle through the 
earth's atmosphere depends largely on l lie pro- 
vision that is made for reducing aerodynamic heat 
transfer to the structure of the vehicle. Analyses 
of tiie heating experienced during reentry have 
been made for both ballistic vehicles (ref. 1) and 
manned vehicles (ref. 2). 

For the purpose of the present report the types 
of reentry are categorized as follows: 

(a) Lifting vehicles of constant lift-drag ratio 
which reenter the atmosphere at very small angles 
(so that skipping does not occur) and which ex- 
perience maximum decelerations less than about 

fy/ 

(b) Xonlifting vehicles which reenter at small 
angles and which experience maximum decelera- 
tions between S// and 14 y 

(c) Ballist ic vehicles which enter at larger angles 
and experience maximum decelerations greater 
than 14// 


The vehicles considered in (a) and (b) are suit- 
abb' for manned reentry, whereas the decelera- 
tions associated with higher entry angles (typo 
(<*)) generally exceed human tolerances. 

The heating experience of the manned vehicles 
also differs from that of the higher-entry -angle 
ballistic vehicle. For tin*, ballistic vehicle the 
maximum heating rates are such that surface 
temperatures may exceed the melting temperature's 
of metals which have been considered for heat- 
sink shields. For manned vehicles the flight 
duration is much longer and, although the maxi- 
mum heating rate is much lower, the total heat 
input exceeds that of tin* ballistic- vehicle and the 
use of a metal heat-sink shield becomes inefficient 
from a weight standpoint. 

As an alternative to the heat sink, consideration 
has been given to the use 1 of ablation materials. 
The term ablation applies when there is a removal 
of material (and an associated removal of heat-) 
caused by aerodynamic heating, and therefore 
embraces melting, sublimation, melting and subse- 
quent. vaporization of the liquid film, burning, 
and depolymerization. 

Several approximate analyses have been made 
of the steady-state shielding effects which result 
from this removal of material. Aerodynamic 
melting has been considered in references il, 4, 
and a, sublimation in reference (>, and simultaneous 
melting and vaporization in references 7 and 8. 
A general treatment of the boundary layer with 
mass addition has been given in reference 9. 
The problem of keeping the vehicle structure' at a 
suitably low temperature cannot be answered by 
investigating a steady-state situation, however; 

1 
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consideration must he given to the problem of 
insulation, and the conduction of heat to tin* 
structure is a transient phenomenon. 

The suitability of a heat shield (whether heat 
sink or ablation material) depends on the weight 
required to keep the structure below a given 
temperature, and a simple quantity of merit such 
as the effective heat rapacity or effective heat of 
ablation gives no indication of the severity of the 
insulation problem. The use of high-ablat ■ion- 
temperature materials such as graphite, for ex- 
ample, which has high thermal conductivity, 
could haul to an intolerable heating condition 
although the effective heat capacity is higher than 
that of most materials. It is possible that the 
high-ablat ion -tern perat lire mat ('rials experience 
less mass loss because of radiation cooling, but 
this is not necessarily desirable since the high 
surface temperature makes tin* insulation problem 
more severe. 

It is seen, therefore, that the problem of main- 
taining a cool vehicle structure is twofold: there 
must be adequate provision of material for abla- 
tion, and there must be sufficient insulation to 
prevent the structure from becoming hot. 

Tin* effectiveness of an ablation material de- 
pends on its capability to dispose of heat by 
convection in the liquid film, as latent heat, and 
by convection in gaseous form in the boundary 
layer. It has been shown (refs, (> and 7) that in 
general an ablation shield is most effective when 
a large fraction of the mass loss undergoes vapor- 
ization. When sublimation takes place there is 
no liquid film; all the mass lost from the shield 
undergoes vaporization and subsequent- convec- 
tion in the high-temperature boundary layer and 
thereby removes a large amount of heat. For 
this reason a material which undergoes sublima- 
tion rather than melting is generally more effi- 
cient (apart from considerations of latent heat). 

Naturally, tin* choice of ablation material will 
be dictated by the type of vehicle and its heating 
history during reentry. For a nonlifting vehicle 
whose dimensions are such that the heating rates 
experienced are too high to lx* balanced bv radi- 
ative cooling, it is advantageous to use a material 
with a low ablation temperature in order to 
reduce l lie insulation problem. For lifting ve- 
hicles, which experience lower heating rates over 
most of the surface, the primary means of cooling 


would be radiative except at the leading edges 
where the limited use of a material with a high 
ablation temperature would seem more appropri- 
ate; ablation would then take place only near the 
peak heating condition. 

The materials here considered most suitable for 
the reentry of a manned nonlifting capsule arc 
therefore those which undergo sublimation at a 
low temperature (say less than 1,500° R) and 
have low conductivity so that no further insula- 
tion is required. The absence of a liquid phase 
insures that the material is removed in gaseous 
form and therefore eonvects a large amount of 
heat from the shield, and also precludes the pos- 
sibility of liquid-film instability. 

The purpose of the report is to develop an 
approximate' method of solution for such a shield 
from which may be determined the* total sub- 
limation of material during reentry and the 
temperature distribution within the remaining 
shield; the ablation temperature is such that 
radiation may lx* neglected. The analysis is 
directed toward obtaining results useful for 
engineering purposes. 


SYMBOLS 


A 

reference area for drag and lift, sq ft 

a 

horizontal acceleration, <j units 


f D 

c 

drag coefficient, 
specific heat, Btu/(ll>) (°R) 


t'n 

specific* heat at constant pressure, 


Btu/(lb)(°R) 


a 

gravitational acceleration, ft /see" 


/ A// 

effective heat- capacity, Btu/lb 


h 

altitude, ft 


J 

mechanical equivalent of heat, ’ 

778 


ft-lb/Btu 


k 

t hermal conductivity, Bt n/ (ft ) (sec)(°R) 

h 

ratio of local heat flux to that 

at 


<7o 

stagnation point , 7 - 


k u 

average* value of beat flux relative* 

to 


stagnation-point value, ~ 

fijs t/xp.O 

dS 

L 

latent heat of sublimation, Btu/lb 


LID 

ratio of lift force to drag force 


! 

characteristic length of vehicle, ft 


M 

mass of vehicle, slugs 


m 

mass ablated per unit area, Ib/sej 

ft 

n Pt 

Fraud tl number 
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Schmidt number 

total convective heat absorbed per 
unit area, Btu/sq ft 
dimensionless heat absorbed per unit 
area 

local convective heat-transfer rate per 
unit area, Btu/(sq ft) (sec) 
dimensionless heat-transfer rate 
radius of curvature of nose, ft 
Reynolds number, p«,V7//i „ 
distance from center of earth to orbit, 
ft 

surface area wetted by boundary layer, 
sq ft 

temperature, °R 
t ime, see 

1/2 r 

27 

tangential velocity component, ft /sec 
circular orbital velocity, (</r) 1/2 = 20,000 
ft /see 


total veloeitv, ft /sec 

. T-T w 

temperature ratio, wv 

-*■ a ~ oo 

weight of vehicle at earth’s surface, lb 
outward normal distance from initial 
position of ablation surface, ft 
dimensionless function of u determined 
by equation (4) 

outward normal distance from ablation 
surface, ft 

fractional temperature rise of gaseous 
material 

atmospheric densitv decav parameter, 
ft- 1 

dimensionless ablation rate 
dimensionless heat content 
fractional decrease in mass loss due 
to latent heat 

integral thickness of heated layer in 
solid shield, ft 

latent heat parameter, equation (IH) 
or (34) 

coefficient of dynamic viscosity, slugs/ 
ft -sec 


dimensionless distance 



1/2 


density, slugs/cu ft 


p* value of reference density for exponen- 

tial approximation to density-alti- 
tude relationship, 0.0027 slug/cu ft 
r dimensionless time, tjt f 

<p flight-path angle relative to local hori- 

zontal direction; negative for de- 
scent, deg 

X conduction parameter, equation (04) 

Subscripts: 

0 no sublimation 

00 free stream 

1 gas produced by sublimation 

2 air behind shock wave 

a sublimation condition 

b solid shield condition 

e external to boundary layer at stagna- 

tion point 

f final conditions 

i initial condition 

surface condition 
xp stagnation point 

Superscripts: 

differentiation with respect to 17 
dimensionless quanlitv 
~ mean value 

ANALYSIS 

ASSUMPTIONS AND APPROXIMATIONS 

The analysis of unsteady sublimation of material 
from the nose of a body during reentry into the 
earth’s atmosphere requires a knowledge of the 
heating experienced by the vehicle. Throughout 
this report the equations which describe the motion 
and heating history of the vehicle are those de- 
veloped in reference 2; the present analysis there- 
fore contains all the assumptions stated in refer- 
ence 2 in addition to those stated now which 
concern the properties of the material and its 
behavior during sublimation: 

(1) The density p bf thermal conductivity k by 
specific heat c by and latent beat of sublimation A 
have constant values 

(2) The Brand tl number N Pr and the Schmidt 
number N Sc of the gas mixture in the boundary 
layer and the specific heats c p l and are constant 

(3) The sublimation temperature T a is constant 

(4) Sublimation leaves a smooth surface and 
causes negligible change in the shape of the shield 
compared with the scale of the noses 

(5) The boundary layer remains laminar for 
most of the reentry 
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To facilitate a simplification of the conduction 
j)rol)lcin the following approximations art* made: 

(0) The material is sufficiently thick to allow 
use of the “infinitely thick slab” approximation 

(7) Conduction of heat takes place along lines 
normal to the surface 

A lew remarks in justification of the assumptions 
and approximations seem relevant. Assumption 

(1) of constant properties may not always he 
appropriate, hut it is hoped that the analysis will 
give useful results when mean values ot p f) , c f >, and 
k h are used. The effect on the gas-layer shielding 
of variations in X f > r and A, Sf is secondary compared 
with variations in and c p » (see assumption 

(2) ); however these specific heats appear only in 
the ratio <* p-1 /c Pi2 , which should have little varia- 
tion with temperature. 

The sublimation temperature T a is correctly 
determined by consideration ot the phase equi- 
librium between the solid material and its vapor; 
in practice, however, T a varies within a limited 
range' only and assumption (H) is justified. The 
assumption of smooth sublimation is actually a 
required property of the material it it is to be a 
successful shield. 

The Reynolds numbers experienced during re- 
entry depend on the parameter H /C y ^l for the 
vehicle. For manned-capsule reentry, typical 
values of the Reynolds number indicate laminar 
flow throughout the flight through tin 1 upper 
atmosphere where the most severe heating is 
experienced. 

MOTION AM) H K ATI Mi DURING RKKNTRY 


veloeit >' romponi'nt u ~ ~ aIu ^ a 1U ‘ W bide- 

pcndci t variable 4 





C 7 .l 


(a) 


in t he following form : 

HZ' ' - (z ' -|)= 1 cos’ v- (fir) ' ~ cos* <p (4) 

For reentry at small angles th.e tangential com- 
poncin of velocity u is approximately equal to 
the Ural velocity, and cos \ . 

The initial conditions used lierein are those 
appropriate to entries starting from a circular or 
near-cir(*ular orbit at high, altitude 4 : 

X, 0 (5a) 

Z/ - (^/‘) ,/J sinv? f (5b) 

at u - l . 

All quant it ies of interest for example, deceler- 
ation, luxating rate, and Reynolds number can 
In 4 det ermineel from the solution of equation (4). 
The deceleration is given by 

— a — -f - (fir) ]r2 UZ - iW TiZ ((>) 
<1 (ft 

and ti e elapsed time by 

t-r, (fry)" 1 ' 2 / = 27 (uZ)~ l (ITi (7) 


An analysis of shallow reentry into the earth’s 
atmosphere, for both lifting and nonlifting vehi- 
cle's, has been made in reference 2. The results 
of inteivst for the present application are* included 
here* for completeness. 

The variation of density with altitude* is 
assumed to la* exponential in form; that is, 



where p* and have* the numerical values 

p* ---0.0027 slug/eu ft (2a ) 


The free-stream Rey nolds number pe*r unit length 
is 


IL 

f 




( 8 ) 


As indicated in reference 2, the Reynolds numbers 
cxjieri *nced during satellite reentry arc* such that 
laminar How may be expected to occur for the 

greater part of the reentry ^typical values 

-5-^40 and Z=0.5 give! — .® — 1 AX 10 ,r ’ V 
( d A ( / 


l 

2df>0() 


ft 


(2b) 


The* equations of motion are combined and 
written in terms of a dimensionless tangential 


The stagnation temperature of the stream is 
given by 
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,vl.('iv 7, 4M2° H. 'Pirns 


The I aminar convective luxating rate which 
would be* e*\perie*need by a nonablating body is 

r \f * - 

given in terms ot >, . , .* a, and Z as 

f i>Ali 




'/n = «*- 


e»AT*— 7’»)"| ry ii4 

1 n* • (w/) 

2 J(j _ 


h\= 1 - is a factor determined bv the shape of the 
nose of the body, and the viscosity law \x t u ^ u 

is used. The term is small compared 

1 v„- 


witli t/“ and is neglected ; thus 

q„=u\uZ)' :i (11) 

The* constant K in equation (10) is that used in 
reference 2; that is, K— 500. Then 


q t) =~M)k 


M .,X : AmZ) lri Bt,u/(sq ft)(swi) (12) 




Tin* total heat which would be accumulated by 
a nonablating shield is found by integration of 
equation (12) to be 

/ \f \i/2 _ 

Qo= lo,900/f n ( ( ' }> \i{J V" n O'*) 


Qu=J”* H 2 {uZ) h ‘du (14) 

and 

f J^'IS (15) 

*■> »/.s <Jsp, n 

The* foregoing relations (eqs. (t>) to (If))) are- used 
in the determination of tin* mass required for 
sublimation and the accumulation of heat within 
tin* solid shield. 


GENERAL EQUATIONS FOR SUBLIMATION AND HEAT 
ACCUMULATION 

Before the vehicle reenters t lie atmosphere the 
ablation shield is assumed to have uniform tem- 
perature 7\„ . In the early part of reentry tin 1 
shield is heated until the surface temperature 
reaches the ablation temperature T„. During 
this preablation heating period tin* problem of 
conduction of heat through the shield can be 
treated without difliculty since 1 the heating rate 
is known (eq. (10)). When sublimation of ma- 
terial occurs, however, the conduction of heat 
within the material depends on tin 1 rati* of mass 
loss from the surface 1 , the problem becomes non- 
linear, and the exact solution involves lengthy 
numerical procedures. (See, for example, ref. 
10 .) 

It is not the purpose of this report to obtain 
exact- solutions of the nonlinear equations; rather, 
approximate* results are obtained which show all 
the important parameters that (‘liter into the 
problem and give estimates of the material re- 
quired for sublimation and for absorbing the 
heat conducted to tin* interior. Figure* 1 is a 
diagram of the* heat shield under consideration. 

The development- of the* differential (‘([nation 
and appropriate* boundary conditions is given in 
the appendix. For the* present it is more* useful 
to eonsideT a heat energy balance* as follows: 

Q(t) = fc,(7;-7h)f /dm + P^f 0 (7- 'O/s 

(lb) 

Net beat Heal absorbed by Heat accumulated by 

input at sublimated material remaining material 

surface 



Figure I. Sublimation at the nose. 
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<) 

When an integral thickne'ss 6 defined as 

*(/) = ] r/5 (17) 

is introduced! into equation (Hi) the following 
equation results: 

Q( t ) -- [ (‘b( l a — / 0 ,)+/d M ~b P //*/»( I* I « ) 0 (IS) 

Alternatively, (‘({nation (IS) may he written) in 
differential form as 

q(t) \ e„(T a - T,)- \ I A +p*r. j f [('/’,- TJ0\ 

(19) 

An additional ('({nation (a boundary condition 
at the surface' of sublimation) is written 

' li,] ~ L )//' + ( A> dO;-0 

(20) 

Net heat-iransfVr Rate of heat al>- Heat-transfer rate 

rati* to surface sorbed in phase to interior 

change 

It is important to note that in equations (19) 
and (20) the lieat-transfer rate q(t) is that which 
tilt 1 shield actually experiences and is itself a 
function of the rate' of sublimation. Throughout 
this report the quasi-steady relation for the 
reduction in heat-transfer rate dm' to the intro- 
duction of mass into the 1 boundary layer is used 
(se'(‘ ref. 0): 

(2i) 

wlmre* #/,)(/) is the heat -transfer rate experienced 
hv a nonablating hodv at the surface* temperature 
Ta- 
in equation (21) c p is the effective mean specific 
he'at and a(T c —T u ) the effective temperature rise 
of the mass converted in the boundary layer. 
The* expressions for a and c p derived in reference 
6 for a laminar boundary layer are* 

«=1_‘ A', (22) 

and 

— ( 23 ) 


where* v? is the effective concentration of the 
shield mat e*rial in gaseous form in the* boundary 
laye'f and is given as a function of N Sc in reference* 
f>. 

The unknown luxating rate q(i) is eliminated 
from eej nations (19) and (20) by use* of ('({nation 
(21) to give 

— [( h(T u — T m ) -\-L-\-ac p ( 7 t: — I „)] jy 

+P^ f l(T-Tj6} (24) 

and 

< ln u)-\L \ac„(T- t u )\ w 


A cor iparison of (‘({nations (24) and (25) shows 
t hat 


( 


h 


dT\ 


c»(T a -rj 


dm 

df 


+ phCh -jj [ ( — T , ) 6] ( 2 () ) 


which, vlum integrated, give's 




d! Ch( T a - T M )m -\ p h c„(7\— T^)d 


(27) 


Equation (27) shows that only part of the heat 
transferred from the surface toward the interior 
is arm nulated (an amount p*<u,(7\.— 7V10), t he 
remainder be*ing re<{uired to raise* the* mass rn 
to the* r hlation temperature. 

In oided 1 to solve* equations (24) and (25) the 


relatioi ship he*t .ween 



and 6 must 


be* 


known the* exact, determination of this relation- 
ship would involve the* solution of the nonlinear 
coneluc ion equation (as shown in the appendix). 
In this report only approximate relations between 



and 0 are us<*d. 


Befo e* attempting to take account of the con- 
duct ior of heat within the* shield it is useful to 
make a simple* analysis of the sublimation to obtain 
e*xpre\ss ons for the* total mass loss during reentry. 


SUBLIMATION OF MATERIAL FROM THE SHIELD 

The purpose* of an ablation shield is to reduce 
the he* it-transfer rate at tin* surface from the 
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aerodynamic* rate c/ (( to a value ^ ) by pro- 
viding material \vbic*li absorbs heat (through latoni 
heat of sublimation) and converts heat in the gas 
boundary layer; this situation is refleeted in (‘((na- 
tion (25). If this process is successful, then 

(jc b ) <C (jo for most of the reentry and an 

upper limit to the rate of mass loss can be ob- 
tained by neglecting heat, conduction in the solid. 

From equation (25), with yk b ^ - ) =0, 


d fn_ q t) 

(If L + ar ; ,( /).— /),) 


(2S) 


Alternatively, if it 
accumulat ion of heat 


assumed that the rate of 


£ 

(It 


\p h c h {l\—T^)d\ 


is small compared with the rate of disposal of heat. 

k„( 7*. ) — /,-!■ ac „('!]. - T„)\ ™ 
equation (24) becomes 


(I m qo _ 

df o,(7 n — / ut! ) + 7v-(~ or p( 7 c — 7 a ) 


(29) 


Equation (29) is the ((uasi-steadv (expression for 
the rate of mass loss but is not , in general, an upper 
limit. Substitution for r/ () and T t — T a ~ T e — r i\ 
from (‘((nations (12) and (9), respectively, and use 
of the relation 


(I m dm du 
dt du dt 


1 _ f . d m 


= - ^ uZ 


du 


m 


( u f 

du c„ du 1 


00 


whei 


and 


S~0-4x) 




L_ 

1 u/ V 

2 J» “ <v 


m 


(33) 


when (‘([nation (28) is used, or 


X- 


L+r f) (T (l -TJ 

1 V <‘r 

2J9 a e,.* 


Maximum internal shielding enthalpy 
Maximum (external shielding enthalpy 


(84) 


when (‘([nation (29) is used. 

Equation (81) shows immediately tin* impor- 
tance of the parameter 



Heating coefficient 
Gas shielding coefficient 


which depends on the vehicle size and shape 
/ M Y /2 

through y and the properties of the gas 

boundary layer. It is scam that, the mass loss can 
be reduced by ({(‘signing the vehicle so that 
M/(' i} Ali is small (low mass, high drag, blunt nose) 
and by choosing an ablation material having a high 

shielding coefficient a ( v ■ The factor a depends 

G ,2 

on the nature of the boundary layer and is given 
in equation (22) for a laminar boundary Inver. 
The specific-heat- ratio (‘ i Jc l>:2 is given in equation 
(28) in terms of c^jc,^ and it is seen that this 
ratio should be large in order to reduce 1 the mass 
loss. 

A second parameter, tin* enthalpy ratio 

^ /> Shielding due to latent heat 

1 u / e» Gas-lav( k r convective shielding 

shows the effect, of the latent heat L in reducing 

the mass loss; when a ( !> is small the (‘fleet of 

G.a 

1 a t ('ii t heat becomes more important. 

Mass loss. — An upper limit to the total mass 
loss is obtained by integration of equation (81) as 


where 


m— 1 .18 b 


6 


M V 



m lb/sq ft 


w 




(:w) 


(36) 


r.t::r.ni t;i 
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and Ti (l and 77/ arc 1 , respectively, tin* values ol a 
when sublimation begins and ends. 

Kquation (25) with J* n = 1 gives the mass loss at 

the stagnation point, and the factor A'n— (IS 

modifies this mass loss according to the variation 
of heating rate over t lie surface of the shield. The 
analysis of reference 2 does not apply at the condi- 
tion 77=1 and it is necessary to assume a value 1 
77<T 1 as the upper limit of 77. The nominal value 1 
— 0.995, used in reference 2, is also used here 1 as 
the value 1 at which sublimation begins. The 1 
lower limit depends on the 1 ablation temperature 
of the material, since 1 ablation will cease 1 before the 1 
stagnation temperature of the 1 stream falls below 
the 1 ablation temperature. The value 77/= 0.05 
(which corresponds to a stream temperat ure 1 e>l 
about 200° F) is useal throughout this report; this 
value 1 is considered sudieient ly low to include 1 any 
material now under consideration. 

Thus 777 can be 1 wri 1 1 en 



= n.\ = ,( 1 — t?) 

07) 

where 

y»0 10J5 



“x.„- (TiZ)- l/2 ,lli 

m 


J 0.05 


and 

X 



- . (TiZ)~ 1 " J dii 

>r~\ 

(:s<)) 


(tfZ) _1/ V« 

J 0.05 


Kepiation 

(28) shows that, even wlien X 

0 (do 


mass loss, wliereas equation (.29) give's 7?(\), the 
fractional decrease 1 in mass loss due to latent heat. 
Tlie evaluation of the integrals in equations (28) 
and (29) ran be 1 carried out when the 1 appropriate 
Z-functiems are 1 inserted. Tlie 1 dependence of 
these 1 functions on FJl) and — <p t has been discussed 
in detail in reference 2. 

Effective heat capacity. -When comparing an 
ablation sbielel with a soliel “heat-sink" shield (for 
exxample, copper eu* beryllium) it is convenient to 
int rod uce an effective heat capacity defineel here by 
the 1 following ratio: 


ir 


Total heat which would be 1 absorbed 
bv a nonablating shield 


Total mass loss from ablation shield 


The total heat absorbed Q a is given by equation 


(12) as 


( K \i— 15/JOOA‘n 


(j'uAji) fl 


when 1 , from equation (14), 


_ /*<>. 995 

J 0.05 


u 2 (uZ)~ ir2 du 


and 1 he^ mass loss is given by equations (25) 
and (2b). 


When 
result is 

these expressions are used the following 
obtained: 


//.„= 1 3,500a ' " ll e „ Btu/lb 

Cp , 2 

(40) 

win 1 re 1 



and 

^7.//=(77,.//) x _ 0 (i— ^) _l 

(4 In) 


/'o.y05 

I 7r (77Z) ~ ]J2 </u 

(4U>) 

{TtZ)~ m dTi 

J 0.05 

Kqualie ns (40) and (41) show that // f// e lepends 
only or the 1 properties e >f the ablation material 
(throng i c f Jc Pt 2 and X) and on the vehicle trajectory 
(since 1 (1 is a function only of trajectory). 

Relations between mass loss, deceleration, and 
time of reentry -Jt is seen from equations (20) 
and (22) that the 1 maximum sublimation rate 
depend}- on the value of X, since 1 


(I m 
dt 


a 2 

P+X 


(TiZ) m 


For large values of X this maximum occurs when 
U 2 (uZy 2 is greatest that is, at peuk heating 
(see eq (11)) and when X---0, it occurs when 
(77Z) 1/2 s greatest that is, at peak deceleration 
(as see i from eq. (b)). Thus, in general, the 
maximi m sublimation rate occurs between peak 
heating and peak deceleration. 

In general, the total mass loss will depend on 
the tot d lime taken to complete reentry since 
777 x= ,i is a function of the 1 trajectory. It is of 
extreme interest, therefore 1 , to determine how tlie 
mass lens may be reduced by allowing the vehicle 
to complete reemtry in a short period of time but 
with the reservation that the maximum decelera- 
tion be kept to a tolerable level. 
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The relationship between the total mass loss, 
horizontal deceleration, and time of reentry 
becomes apparent when equations (<>), (7), and 
(38) are recalled: 


and 


—a= (pr) V2 uZ 



Since uZ^{uZ) max , it is seen that 


Mx=o^ 



— 1/2 


and since 




( uZyhtv 


(by a simple application of the Schwartz integral 
inequality), then 


Tims, in general, /Mx-n satisfies 



(where the numerical values of (j 8/‘) ly2 and (j tty) 1 * 2 
have beam inserted), a result independent of 
vehicle characteristics or trajectory. When an 

average value of uZ equal to 9 (JlZ) ma . r is used, the 

following simple rule may be expected to hold: 


m\=i i 



The foregoing relations between mass loss, 
deceleration, and reentry time provide very simple 
expressions for the sublimation during reentry. 

Application to particular vehicles. -The simple 
approximate rule A is investigated by comparing 
exact results for w x=0 (both analytic and numer- 
ical) with equation (42). The integrals that 
appear in equations (38), (30), and (41b) are 
evaluated for the following types of vehicle: 

(a) Lifting vehicle at zero reentry angle 

(b) Nonlifting capsule at small entry angle 
(e) Ballistic vehicle 

Although the decelerations experienced by type 

(c) vehicles usually exceed human tolerances, it is 
nevertheless of interest to compare the results 
with those for types (a) and (b). 


(a) Lifting vehicle ^>1;— #>*= 0 ^ 


The appropriate Z-function is written 
1 — u 2 


Z- 


50 f } v 


143) 


From equation (38), 

( r\ \n /'o.ws 

m v) L 

1.418 (:i() (44) 

and from equation (39), 


*1= 



(l — ?/ 2 ) 1/2 </u 
-u 2 )- ]/2 (lu 


(,-x) («> 


That is, 

Kale A: The total mass loss ran es inversely as 
the square root of the maximum horizontal decelera- 
tion. 

Alternatively, comparing the integral expres- 
sions for t and 7a x =o leads to the following relation: 

mx =0 ~ ( 2 V) 

That is, 

Rule B: The total mass loss varies as the square 
root of time for reentry . 


^ When the limits of integration are 0 and 1 , 
/ X V /2 \ 

77 = L ^ j ■ j From equation (41b), 

0.995 __ _ 

«•'« « (4(») 

(i -u 2 y in <m ~ 

Jo. 115 

(Again, wlion the limits art* 0 ami 1 .) 
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The maximum deceleration is 

— •==-:*I)( : mZ)„ i , — (^) (47) 

am{ e<{uation (42) gives the approximate result 

Wx.o=2'«(»()^y (48) 

A comparison with equation (44) shows that 
the error is loss than one-half of 1 percent. 

00 Nonlifting capsule ^=0;— <^<5°) 


and equation (42) gives 


Wx=o — 1 8.80[80 sin( — ^?,)] 2 (52) 

A comparison of equations (50) and (52) shows 
that 1 1 o * error incurred by the use of the approxi- 
mate relation (42) is less than 1 percent. 

Evaluation of the integrals in equation (41b) 
gives 


iJLff) 


0 



(58) 


I he Z-funetions tor nonlifting vehicles at small 
reentry angle 0^ — ^f\5° are tabulated in refer- 
ence 1 2, together with associated functions of Z. 
These functions, which are included in table I of 
the present paper for easy reference, are used to 
evaluate nt x=n% ??, and (IJ tff ) K=n by numerical 
integration and tlie results are given in table II, 
together with those tor cases (a) and (e). For 
case (b) the agreement between the results for 
hy rule A and the exact values is again 
l 0 

good except lor — , where the error is 10 
percent . 

Tlit* <|ii!i n t it it*s (7?, // ) x=l , and i? for ease (l>) agree 
almost exactly at — with th<> corresponding 
values for the lifting vehicle. As — Vi increases 
from zero, however, both (7? f// ) x=n and v change, 
hut approach constant values independent of <p, 
when — <pi>o° as shown by t lie following eon- 
shlerat ion of ease (e). 


which is independent of — <p u 

The junction 77(A) is also independent of — <p t 
lor <pi^> 5° and is evaluat-e<l numerically from 
equation (89). Table II shows that yj-> 1 as 
X ->oo for all east's. Also, as X — >00 t t he effective 
heat rapacity ll iff must tend to a limit A, the 
latent Imat, or 

II cff =(IT t . rr ) x ^( l-,)- 1 --* 

This behavior may be verified for east' (a) where 

<7b,W. -,)-,([, 

It is seen that because the exponent of 

I \ ‘qtnd to IL iff ^ This behavior sug- 

gests that t?, which is equal It) for 

case (a), should be more generallv 


(c) Ballistic vehicle ^^=0;— v?,>5 0 ^ 

The retpiired Z-function is now 

Z=80 sin(— ^/)?7(— log,. u) (49) 

which is inserted into equation (88) to give 

n 0.995 

-o=[3() sin(— v?,-)] {—u 2 log r u) m dn 

J o.or> 

=».»2 [30 sin(-^)] m 




(50) 

The maximum decelemtion is also easily obtained 
as 


_(_Ulu\ _30 2 sin (—#>,) 

V 7 ) <lt /mar 2e 


(51) 


/ X V"r tr \.. a 

anti in fact a comparison shows that this result 
agrees w th those in table II to within 1 percent. 

The foregoing results for 7 (77,//) x= , } , and 
V are shown in figures 2, 8, 4, and 5; in figure 2 
the values of t he maximum deceleration required 
for equation (42) were obtained from figure fi(a) 
t)f re fere ice 2 for 0^ — <£,<^fi°, and from equation 
(51) of t ie present report for —<p>>b 0 . 

In fig ire 8, the values of the maximum de- 
celeratio i were taken from figure 10 of reference 

2 for and from equation (47) of the 

present report for ~>1- 
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TAKLK 1. -VALIJKS OK Y AM) RKLATKI) QUAXTITIKS FOR /,//>= 0 AND 



7 

qu 

(Jo 

/, sec 

7 

5o 

Qn 

/, sec 

7 


Qu 

l, sec 

u 


— tp r 

0° 





1° 

2 



- 

1° 


0. 995 

0. 00058 

0. 0288 

0 


0 

0. 00181 

0. 0857 

0 

0 

0. 0020 

0. 0504 

0 

0 

. 99 

. 00105 

. 0890 


150 

140 

. 00270 

. 0507 

.118 

72 

. 0058 

. 0710 

. 080 

87 

. 98 

. 00407 

. 0050 


888 

288 

. 00008 

. 0788 

. 205 

140 

. 0108 

. 0988 

.191 

78 

. 90 

. 01815 

. 1055 


550 

809 

. 01457 

. 1090 

. 400 

199 

. 0229 

. 1800 

. 842 

108 

. 94 

. 0241 

. 1880 


088 

841 

. 0258 

. 1808 

. 598 

229 

. 0801 

. 1028 

. 451 

128 

. 92 

. 0809 

. 1559 


792 

800 

. 0878 

. 1578 

. 095 

248 

. 0505 

. 1824 

. 587 

111 

. 90 

. 0515 

. 1744 


875 

874 

. 0519 

. 1751 

. 777 

202 

. 0000 

. 1974 

. 009 

152 

. So 

. 0989 

. 204 l 

1. 

080 

890 

. 0984 

. 2080 

. 982 

284 

. 1092 

. 2201 

. 750 

170 

. 80 

. 1485 

. 2108 

1. 

140 

411 

. 1421 

. 2158 

1. 042 

298 

. 1580 

. 2275 

. 858 

188 

. 75 

. 1991 

. 2174 

1. 

228 

421 

. 1970 

. 2 1 02 

1. 120 

809 

. 212 

. 2248 

. 988 

192 

. 70 

. 200 

. 2090 

1. 

288 

429 

. 257 

. 2078 

1. 191 

817 

. 270 

. 2180 

. 990 

200 

. 05 

. 524 

. 1959 

1. 

840 

480 

. 82 1 

. 1980 

1. 248 

824 

. 888 

. 1900 

l. 047 

207 

. 00 

. 892 

. 1740 

1. 

881 

442 

. 889 

. 1740 

1. 285 

880 

. 898 

. 1759 

1. 088 

2 1 8 

. 55 

. 408 

. 1520 

1. 

415 

448 

. 400 

. 1521 

l. 818 

880 

. 400 

. 1581 

1. 122 

218 

. 50 

. 580 

. 1294 


442 

458 

. 588 

. 1291 

1. 840 

841 

. 587 

. 1295 

1. 149 

228 

. 45 

. 010 

. 1001 

1. 

404 

458 

. 007 

. 1058 

1. 807 

840 

. 008 
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1. 170 

228 

. 40 

. 084 

. 0887 

1, 

481 

408 

. 081 

. 0885 

l. 885 

85 1 

. 080 

. 0884 

l. 188 

288 

: . 55 

. 757 

. 0081 

1. 

495 

408 

. 754 

. 0029 

1. 899 

850 

. 752 

. 0028 

1. 202 

288 

. 80 

. 827 

. 0448 

1. 

505 

478 

. 824 

. 0447 

1. 409 

801 

. 821 

. 0447 

l. 212 

244 

. 25 

. 892 

. 0295 

1. 

518 

479 

. 890 

. 0295 

1. 417 

807 

. 885 

. 0294 

l. 220 

249 

. 20 

. 949 

. 0174 

1. 

519 

480 

. 947 

.0174 

1. 428 

874 

. 942 

. 0174 

1 . 220 

250 

. 15 

. 992 

. 0087 

1. 

528 

498 

. 991 

. 0087 

1. 427 

882 

. 980 

. 0080 

1. 280 

204 

. 10 

l. 009 

. 0082 

1. 

525 

504 

l. 009 

. 0082 

1. 429 

892 

1. 005 

. 0082 

1 . 282 

275 

. 05 

. 958 

. 0005 

1. 

527 

521 

. 958 

. 0005 

1. 481 

409 

. 950 

. 0005 

1. 284 

292 

| . 025 

. 825 

. 0001 

1. 

527 

587 

. 825 

. 0001 

1. 481 

425 

! . 824 

. 0001 

1. 284 

807 

j 


— *fi 

2° 




— <p, — 

8° 


! 

— <pi- 

•r 


0. 995 

0. 0052 

0. 0712 

0 


0 

0. 0078 

0. 0872 

0 

0 

0. 0105 

0. 1012 

0 

0 

. 99 

. 0105 

. 0999 


057 

18 

. 0157 

. 1222 

. 047 

12 

. 0209 

. 1410 

. 040 

9 

. 98 

. 0210 

. 1878 


180 

87 

. 0818 

. 1082 

.111 

24 

. 0410 

. 1989 

. 097 

18 

. 90 

. 0422 

. 1855 


240 

55 

. 0024 

. 2250 

. 201 

87 

. 0827 

! . 2597 

. 175 

28 

f . 94 

. 0088 

. 2104 


827 

00 

. 0984 

. 2018 

. 208 

44 

. 1 288 

. 8008 

. 288 

88 

. 92 

. 0858 

. 2878 


898 

74 

. 1248 

. 2802 

. 828 

50 

. 1088 

. 8280 

. 281 

88 

. 90 

. 1080 

. 2525 


449 

80 

. 1550 

. 8025 

' . 870 

54 

. 208 

. 8402 ! 

. 822 

41 

. 85 

. 1 05 1 

. 2700 


500 

92 

. 28 1 

. 8201 

. 408 

02 

. 800 

. 8048 

. 404 

47 i 

. 80 

. 224 

. 2709 


040 

100 

. 800 

. 8100 

. 580 

i 08 

. 898 

. 8588 

. 409 

52 , 

. 75 

. 284 ! 

. 2590 


714 

107 

. 880 

. 8008 

. 595 

74 

. 482 

. 4882 

. 521 

50 

. 70 

. 840 

.2411 


709 

118 

. 452 

1 . 2750 

. 048 

78 ■ 

. 508 

. 8089 

. 504 

59 ! 

. 05 

. 409 

. 2178 


815 

119 

. 528 

! . 2404 

. 088 

82 

. 049 

. 2744 

. 599 

08 

. 00 

. 471 

. 1920 


852 

124 

. 591 

! . 2144 

. 717 

80 

. 720 

. 2870 

. 080 

00 

. 55 

. 588 

. 1045 ! 


888 

128 

. 057 

. 1818 

. 744 i 

90 

. 797 

. 2008 

. 055 

09 i 

i . 50 

. 008 

. 1878 1 


908 

i 188 

. 721 

. 1501 

. 707 i 

94 

. 808 

. 1042 

. 070 

72 i 

1 - 45 

. 008 

. li 10 


929 

187 

. 781 

. 1200 

. 787 

98 

. 922 

. 1804 

. 098 

75 i 

i . 40 

. 782 

. 0800 ! 


940 

142 

. 888 

. 0920 

. 802 

102 

. 97 1 

. 0999 

. 708 

79 

. 85 

. 794 

. 0040 ! 


959 

147 

. 890 

. 0084 

. 815 

100 

1 . 018 

. 0781 

. 719 

82 

. 50 

. 858 

. 0455 | 

. 

909 

152 

. 980 

. 0477 

. 824 

1 10 

1. 052 

. 0500 

. 729 

80 

. 25 

. 907 

. 0298 i 


977 

! 157 

. 975 

. 0809 

. 882 

115 

1. 070 

. 0824 

. 780 

91 ' 

. 20 

. 954 

. 0175 


988 

104 

1. 005 

. 0179 

. 888 

121 | 

1 . 080 

. 0180 

. 741 i 

97 

. 15 

. 989 

. 0087 


987 

172 

1. 021 

. 0088 

. 841 

129 i 

1. 079 

. 0091 

. 745 

104 

. 10 

1.002 

. 0082 


989 

182 

t. 010 

. 0082 

. 844 

189 

1. 048 

. 0082 

. 747 

114 

. 05 

. 952 

. 0005 


991 

200 

. 952 

. 0005 

. 845 

1 57 

. 900 

. 0005 

. 748 

181 ; 

. 025 

. 822 

. 0001 


991 

215 

. 820 

. 0001 

. 845 

172 

. 821 

. 0001 

. 749 

140 : 
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TABLE II. TABI LATE!) VALDES OF THE MASS-LOSS PARAMETERS 


Diinrnsimih'ss mass less, m x 0 


Vehicle 
(a) Liflinu, 

Kmry 
; iilT\ 

( 1 :*L’ 

1) 

Kxa t 

, . . /..v L \ 1 J 

Rule A, 

iu-r)r 

, .../ LA 1 '- 

Rule It, 

O.S» ( 7) 1 - 
/ L \ 

(//,//] n 

X = 0 

X = 0.05 

X = 0.1 

X = 0.25 

A = 0.50 

x= 1 .0 

X — CO 





1,4 ( :! %) 

0. 0 

*1) 

*0. 21* 

*0.302 

*0.447 

*0. 577 

*0. 707 

*1 


0 

3.0 4 

2. 71) 

3.00 

0. 500 


u. 22a 

0. ao. r > 

0. 4. r 0 

0. 5*0 

I). 707 

1 


-5 

2. 05 

2 72 

a. 27 

. 4K4 


2;t4 

.am , 

. 403 ; 

. 5*0 

.715 

1 

(ft) N'oulifiinir 

1 

2. 73 

2. 7l) 

2. 70 

. 453 


.24!) 

.337 ; 

. 4*2 

. OOK 

. 730 

I 

< apsulc, 
L 

n = l> 

2 

2. HI 

2. 57 

2. 20 

.412 

<1 

.27a 

. 304 

.511 

. 033 

. 750 

1 1 

a 

2. 17 

2. :«) 

2. oa 

.300 

0 

.2*0 1 

. 3*3 

. 530 

. 050 : 

. 702 

1 


\ 

1.00 

2. ia 

] . M *> 

.ani 

0 

an2 

. 31*7 

. 542 

.000 

. 770 

1 

<e) Ballistic. 

1 ) 


3.32[30 sin{ — sr I]’ 1 -' 2 

a.aojao siuf — vO] -,/2 

a.22[au sin ( — ) |— t/s 

0. 351 

0 

0 . 1:0 

0.43! 

0.572 : 

0. 0*0 ; 

(I. 7*5 j 

1 


v~ 



for case (a); 17 


f x 

O+x/ 


111 ore irciit' rally. 



Fici rk 2. — Dimensionless total mass loss fc r nonlifling vehicles. 
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.5 


1.0 


1.5 

Lift-drag ratio, 


L_ 

D 


2.0 


2.5 


Fici-he 3. -Dimensionless mass loss for lifting vehicles. 
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Initial reentry angle, t deg 

Fkji'kk 4. -Dimensionless effective heat capacity for 
nonlifting vehicles. 


1.04 

X = CO 



u- 

p~ 





v> 8 
v> 

1.00 ’ 

— 

— ( 

!> : 

tn i 

to 

O 

< 

>- — < 

> -o- 

f .6 ( 

a 

o 

> — < 

. j -°\ 

u 

3 t 

TJ „ 

0) 4 

0- — < 

> ■ °1 

o 

O L 


\ 

— c 

j> \ f> - 

t> J 

— 0 ^ 



L ^ i. - n 

S .2 ■ 

Ll 

^ 1 

o 

V 1 rj v 
= 0 , 0 < - <p f < 5 ° 


o ; 

□ 

j 

= 0 ; - <p. > 5 ° 


O 12 3 4 5 6 

Initial reentry angle, - <p f , deg 

Flocre 5.- Fractional reduction in mass loss due to 
latent heat, for lifting and nonlifting vehicles. 


The value's of (//,//)x=o mid V foutui numerically 
for the ease' of jy~() and — — 0 an* in good 

agreement with the* analytical results for yy> 1 
and — (p t — 0: 

W.A.,-5 


ur 


as stHMi in table II, and it is inferred therefore 
that these values are valid for all cases whore 

y } A) and (pi— 0. 

in general it is seen from figure 4 that 0.35 
^(77 ,/A=o ^0.5 for nonlifting vehicles, the low- 
er limit corresponding to — <p,>5°. It is to he 
expected t hat (// (/ /) x=n decreases as — <p t increases, 
since a greater part of the reentry is spent in the 


( 55 ) 
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in table II and shown in figure 0. 

In the foregoing analysis, expressions for the 
total mass loss and the effective heal capacity of 
the shield have been obtained by assuming that 
sublimation starts early during reentry and that 
the accumulation of heat is negligible compared 
with the disposal of heat. Both assumptions lead 
to conservative values for the total mass loss 
(that is, values that are too large). 

An analysis of the heat-conduction problem 
within the shield is desirable in order to justify 
these assumptions and also to estimate the amount 
of insulation required to keep the st nurture 1 cool. 

ACCUMULATION OF HEAT WITHIN THE SHIELD 

The effectiveness of an ablation shield in 


reducing heat transfer to the vehicle structure is 
measured finally in terms of the mass required to 
keep the structure below a given temperature. 
It has been shown in the previous section that 
t he mass loss due to ablation is virtually inde- 
pendent- of surface temperature when 7’,» 1\. 
The mass requirements for insulation, however, 
depend critically on tin* ablation temperature, 
and it is to be expected that tin* use of materials 
with low ablation temperatures will reduce 
considerably the insulation problem with rel- 
atively little increase in the total mass loss. 

In order to estimate the insulation required, 

the unknown relation between (^> and 6 




Eigcue (>. — Variation of mass loss with reentry time- 
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discussed earlier in the 1 report is assumed to be 


(<FLr ,, 'V M 


U'lVTJ 


Although this equation appears at first sight 
to be at most a mule approximation, one should 
expect a relation of the form 


respectively, 

<Zn(/) = [o>(7 T „— 7\J-f L-\-ac — T a )\ 


dm 

~dT 


d 


+ P/A f j-j [ ( 7 x 7 ^ )d\ (o 7 ) 


and 


qM) = {L+«?,JT r --T,,)\ j 7 °°- ( r > 8 ) 
which may be solved for the unknown quantities 


T s and 0 when 


r/ /// 

77 /"" 


-0 (before ablation) or for m 


and 6 when 1\—T a (during ablation). 

Preablation heating. -Before appreciable abla- 
tion occurs (when the mass-loss rate has negligible 
effect on the heat-transfer rate), equations (57) 
and (58) reduce to 


d | f riv r I y \ a 1 ^ />( 7 x 7 , x . ) 

ffo - P/A ^l(/,-/ w )0| = -Q 

and elimination of 6 gives 

(A>Qu=PifJct>(T s — T a y 

When sublimation begins, u U fl and T s 
equation (00) then becomes 

Vu T m ) 2 

which is written in dimensionless form as 

1 . 0 / phCbkt, 


( 00 ) 

= T a ; 

(01) 


«o(M-a)Vo(-W 0 ) = 


kiku AZ__ 

CiM 


(Tu~ rj 2 (10- 7 ) (62) 


(Strictly speaking, sublimation occurs at all 
values of the surface temperature and T s is 


related o dm hit through the phase relation which 
describes the equilibrium of the solid material 
with its vapor. In practice, however, dmjdt is 
negligible 4 except when 7* 1 i os within a limited 
range which includes the mean value T n used 
here.) 

* 1, it is se(*n from (‘({nation ((>2) 


Since 


ku 


(* 

d .r ) 

1 =A(t) h(T - - (5<i) 

,*<) " 

that. 

where A(t) 

is a 

slowly varying dimensionless 

where 

parameter. 

It is 

shown in the appendix that 


N i w 
IIA 

A 




Equations 

(24) 

and (25) are now written, 



r/oO/JVohwJ ~X‘( 1 (r 

(ptfjc h ) h - 


x — - - - 






(«») 

((54) 


Equation (03) determines whether sublimation 
will occur during reentry. If 


«*<X 2 (16 _7 ) 


(<h r >) 


sublimation will not occur, since equation (03) 
cannot be satisfied and the maximum value of 
T s will remain below T a throughout reentry. In 
terms ( f u and Z, the condition that sublimation 
will not occur is written 


X* >10 


Tr(uZ) ]/ - I u\uZ)~ lri du (66) 
J <t J max 


and it \ .ill lx* shown that t lie bracketed expression 
in this inequality is virtually constant (approxi- 
mately 0.2) and is independent of Ljl) and <p t . 
Wild x i s written 


p^AT-T.M h )' 

\Pf)( b/ 


G".*> 


it is sei n that sublimation cannot occur if 

M W 


(l) the parameter 




which deter- 


mines tlie level of the lieating rate, is too small 

(2) the thermal capacity p h c b (T a — T^) is too 
large, or 

/ k b \ ]/2 

(3) the thermal diffusivitv ( — - ) is too large 

v 1 \Ptfib/ 

The materials under (consideration in this report 
have low ablation temperature and low thermal 
conductivity; more specifically, the materials 
under serious consideration have properties with 



AN ANALYSIS OF ABLATION-SHIELD REQUIREMENTS FOR MANNED REENTRY VEHICLES 17 


the following orders of magnitude: p&M)(10 2 ) 
Ib/eii ft. o,«0(l) Btu/(1I>)(°R), /’ h s «0(l(r 5 ) Btu/ 
(ft.) (see) (°lt ), and T tl - r I\ «0(10 2 , 1 0 s ) °K ; for 

vehicles considered here f r~rn ^>10 _1 * Thus 

Oflil/i 

X<Cld a and sublimation will occur. 

It is of interest to determine what fraction of 
t he total flight time passes before sublimat ion takes 
place. The ratio t a jt f is found as follows: 

If t— 0 when 77 — 0.995, then, since 7j 0 is an 
increasing fuii(*tion of 17 during the early part of 
reentry, 

tfo('Wff)>tfn(0.995) 

and 

«oCS a )>7 tt 5o(0.995) 

Therefore 

7»('Wn) Qi)( u «)>7o“(0. 995)7, 

Also 





so that 


(M V n) 

tr 7/" MO-995) 




0)7) 


where the right-hand side of equation (07) i* 
approximately equal to 


Z(0.99f>) 


/Ml. 9!>"t ~|-l 

.995) I (uZ)-'<lu x-(10“ 7 ) 

J0.05 _ 


The ratio /,*.// f is evaluated in a later section by 
inserting the appropriate Z-functions, and it is 

show n that y 1 . 

Total accumulation of heat. In the analysis 
that follows an upper limit- is found for the aceuni- 
mulation of lieat during reentry by assuming that 
the surface of the shield is raised instantaneously 
to the temperature T a at U0. (Although this 
assumption implies an infinite heating rate at 
/-=(), it is nevertheless a, good approximation 
except for small values of t/tf.) 

Equations (26) and (50) arc first combined to 
give 

p '' e m {p, ‘ ey - =P ^~ p ^w m 

and the last term is t lien neglected to give 


Integration then gives the following upper bound 
for 9 : 

P, 0 <(l P f’’ / )' ' (70) 

Tin* accumulation of heat at time t is therefore 


(71) 

and, when / —t /y 

Qf< (-P//7//'/i) 1 “( 7 a — /A ) t / " 

This up})er limit is independent of y^\ n and 

forgiven material properties depemls only on the 
total time of reentry. When Q f is expressed as a 
fraction of (A, / there is obtained 


whei 


Qn,f k n Q iK/ 

- 11/2 

— J ‘ " /• 0.995 

'far u 2 (vZ)~ 17 VS 

J ( 1.05 


(72) 


(73) 


and is evaluated by inserting the appropriate 
Z-functions. 

Application to particular vehicles. -The results 
of tin 1 previous two sections are summarized 
briefly before specific application is made: 

The (plant it v 


X-^ 


XL Y /2 


was showui to be an important parameter of tlie 
heat-conduction problem; numerically X'ClO’ for 
the low-temperature low-conductivity materials 
considered herein. 

Sublimation will not occur during reentry if 
x 2 >l() 7 /q 

whnv 



The ratio t a t f satisfies 


Pb® 


d_ 

dt 


{pb®)< 


Pbkb 

Cb 


(69) 


^<io- W 2 

• f 
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\vll(T(“ 


Z((). !!()")) | (mZ)"WwJ 

The ratio (J/iQ, K t siitisfies 


Qf / 


io- :, x/-3 


where 


„ ,U., ,5Zi H 


j: 


u-(uZ) ]r2 tlu. 


materials; in such a cast; the factors pbfib(T a —Ta>) 


( :) 

\Pl> /,/ 


which appeal* in tin* parameter x 

would be replaced by the analogous composite 
quantities.) 

(2) The ratio 


(70) 


i a Preablation heating period 
1 f Total reentry time 

satisfies 


[''<(&) (10 ,v ) X- 

t f 


, . . 11 . 0.40 1 . \ 

wIhti 1 it is assumed that -r — ~ ~ in eq. (/2) 1* 

A*n 2 / 

Tlie numerical values of F lt F\, and F\ depend 
on the ])articular Z-funetions. In each ease, how- 
ever, the combination of functions is such that the 
total power of Z is zero and tlie factors LjD and 
sin ( <r t ) which (Miter into the (‘valuation finally 
cancel, giving results which are independent of 
L I) and as shown below. 

(a) Lifting* vehicle ^ L “>1;— 0^ 

Wit h TiZ ^ *! , /y, l >, and /^are evaluated as 

F\- 0.2, F h4 ? and t\— 1.1, and are independent 
of LID. 

(b) Nonlifting capsule (-^=0; 0^ — 

Numerical (‘valuation of the integrals gives the 
following average values: /v 0.2, F> 5 0, and 
F,~ 1.5. 

(c) Ballistic vehicle (yy"-();—^,^>5 0 ^ 

K valuation of the integrals by using 
Z o() sin(— pi)u(— loge») 
gives l\ 0.2, F 2 — 14, and /v 1.5. 

It is concluded from tlu‘S(‘ numerical vahu‘s that, 
for cases (a), (b), or (c): 

(1) Sublimation will not occur during reentry if 
X'C > 2XI0 ,i . (This condition does not hold neces- 
sarily if the shield is a composite slab of different 


(:V) The ratio 

Q f Heat accumulated by ablation shield 
Q {1 f Heat accumulated by heat sink 

satisfies 

/Z-<(l-5)(10- :, )x 

Vo./ 

The i ssumptions, made in the analysis of sub- 
limation, that sublimation begins soon after the 
initiation of reentry (when 77=0.995) and that the 
accumulation of heat is small compared wit 1 1 the 

disposal of heat are justified, since y <( 0.05 and 

F < (.15 for x < 10-. 

Vo,/ 

Insulation requirements. -The method of the 
preceding section gives an estimate only of the 
amount of heat accumulated by the solid shield 
at the *ompletion of reentry. The temperature 
distribution through the shield is also of interest, 
howeve *, and an approximate analysis from which 
this distribution can bo obtained is desirable. 
The nonlinear differential equation and boundary 
conditu ns are derived in dimensionless form, to- 
gether with the approximate method of solution, 
in the appendix. The results so obtained are 
summaiixcd briefly and discussed below. 

The emperature distribution at the end of re- 
entry is 
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( M T . . 

*,=<>. 1 7*,, . (1-7,) ( c "- ) * (78) 

<‘n \PtJi f, / 


vcrilic'd that f<l mid f * 0 ~ for large values of e, 

*■*/ 

(in practice, e r >2). 

The mass loss rate din/dt will be small through- 
out reeutrv if the latent heat is large, since 


T— T 

'I'he ratio T y? is shown as a function of £ in 

figure 7. (The more general expression for tem- 
perature distribution during reentry is given in 
the appendix.) 

Tiie accumulation of heat is 


()f= -Tt* TJt (7«) 

7T 

where 

' 4 2 + 2«, Crf 2 

The function f is shown in figure 4 8; it may he 

I 0 r 1 1 T™ - 1 - r r - — 1 






T _ _ 

-~\ - 

•ft 
€ r 

V 

£ 

\° 





\V 








:: 





0 .5 1.0 f.5 2.0 2.5 3.0 3.5 4.0 

Dimensionless distance from ablation surface, ( 

]m<;ire 7. Temperature distribution within the solid 
shield. 


T Poo . t = /p^Y 
T a 7 V s \ktjfj 


.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Average rate of mass- loss parameter, 

Fku'ke 8. — Variation of f with e/. 


L + a n (T-T„) 

<‘,Kt 

Thus, for a low-temperature shield (for instance, 
where x<10 2 ) which also lias sufliriently large 
latent heat, most of tlu v aerodynamic heat input 
is lost in boundary -layer shielding and in the phase 
change, whereas the amount of heat conducted 
from the surface to the interior is sufficient only 
to keep the surface at the ablation temperature 

T a . In the limiting cast 4 when ^ — >0 and 7„=7„ 

throughout reentry, the conduction problem is 
described once more by a linear equation which 
has tin 4 solution (as e ->()) 


T- 

Ta-T" 


4 <«■> 


Equation (80) actually furnishes an upper limit 
to tin 4 temperature 4 distribution, since when abla- 
tion takes place the heat, content of the remaining 
shield is reduced. 

An upper limit to the heat content of the 
material is obtained as 

r o 7 — T 

Q = Ptf'bi I ;x.) j y f y t ~ 


wliieli is simply 


1 2 (pl) r bkl>) ] ’ *( Ai 


At the ( 4 iid of reentry, t = t f and 

Qf ~ 7 3 (p//)W 1/l! ( 7 a / ct ) tf ' * 1 ) 

7T 

This ( 4 xaet limit is about 20 percent less than the 
approximate 4 value given in equation (71). 

The behavior of the solution for large 4 value's of 
e f at the emd of iventry is written 

T n — T,„ 

r / m v 2 1 


= exp 0.1 Ik 


V 

CnAfi) 


-(1-7,)',." 
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Equation (S2) shows the effect of the sublimation 

/ M y« 

VYVIAV 


parameters - 


* /» 

a — 

( p. 2 


and 


rj on 


tempera- 


ture distribution within the shield. 

The insulation requirements -that is, the 
amount of material which should remain whim 
sublimation has ceased in order to maintain the 
structure at some design temperature below the 
sublimation tempera-tun 1 will depend primarily 
on the thermal dilfusivity of the material k h ip f ,c f) , 
the total time of reentry t f , and the ablation 
temperature T i} . The required insulation thick- 
ness can be determined approximately from 
equation (77) or figure 7. 

For engineering purposes the curves of figure 7 
are easily approximated by 


lions have been made, but they an* of such a 
nature as to give conservative results, since upper 
limits have beam obtained for tin* total mass loss 
due to sublimation and for the total heat ao 
cumuh led during reentry. 

It Ins been shown that the total mass required 
for sublimation depends primarily on the para- 
maters 


i , u f 

\('i,A/if _ Healing coefficient 
v p Gas-shielding coefficient 


and 




L\c h (T (l TJ 
1 v* v p 

*j<*y,:* 


T - 7 \, 
T n -7\ 


m 


Maximum internal shielding enthalpy 
Maximum external shii*lding enthalpy 


whi*re 6 is chosen in such a way that the heat This definition of X gives the quasi-steady result, 
content is correct; that is, whereas conservative results are obtained when 


f" 2 

Pb(‘h I ( / ' / oc f) 1 “( l a / ) b( e /) 

Evaluation of the integral using tin* profile of 
equation (82) gives 

(s4) 


Therefore the amount of insulation — pz ahead of 
a station having temperature T is, from equations 
(82) and (84), 


2 / pA h \ l j 

cj 


Z/'yfO/) l<>g <’ 


7 k— 7 V 


7 - 7 ’ 


(Ho) 


The inq)ortance of the shield material parameters 
Pitkbfr'b and T a is evident from the foregoing 
expressions. It is seen also that the amount of 
insulation varies as t f Xrl . 

DISCUSSION 

In tin* preceding analysis, the primary objective 
has been to obtain simple useful expressions to 
describe the sublimation of material from, and the 
accumulation of heat by, a low-temperature low- 
conductivity shield suitable for manned reentry. 

For the sake of simplicity several approxima- 


L 

1 u e 2 c 

2Ju a 7 p 


It is evident that t lie correct interpretation of a — 

( 'p. 2 

is imp< rtant in tin* use of these parameters; the 
quantity arises when the convective shielding in 
the hoi ndary layer is considered and is correctly 
interpreted as 

Enthalpy of gases conveeled 
v { , _ in boundary layer 
c Pr > Enthalpy difference across 
boundary layer 

In reference 7 this ratio has been replaced by 

an emp rical factor 0.08 y yf J ’ whi*re \t is the 

molecular weight, whereas in reference (> simple 
approximate expressions for a and c p /<‘ p ■> w(*ri* de- 
rived under tin* assumption of constant v, iA and 
c „ t 2 an< a[>pear as equations (2*2) and (2)1) of the 
present report. In practice, however, c /M and 
c p ,2 are functions of temperature and the average 
values used should be consistent with the foregoing 
physical interpretation. 
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Furthermore, although the quantity a^\ — 

~ A /V _l, fj is an approximate constant value, t ho 

exact theoretical values for a depart from this 
when the interna! shielding enthalpy A + c a ( 7 t „ — 
T „) is small eom])ared with the gas-layer shield- 
ing enthalpy. (See the discussion of ] OM 
page { of ref. (>.) I his situation may occur in 
the (airly part of reentry if the latent heat of the 
material is small enough, and the use of a value' 
of a that is too large would lead to an under- 
estimation of the shield requirements. 

Although flit' simple expressions (22) and ( 23 ) 
are used herein it is understood that an experi- 
mental verification of tin* ratio 


Enthalpy of gases eonveeted in houndary layer 
Enthalpy difference across boundary laver 

for the particular material under consideration is 
desirable. 

The total mass loss during reentry can be written 

(. m ;f 

m = 1 . 1 S/r j ~ nt\.^ 0 (1 — 17 ) lb/sq ft 

« ( e- 

where m x^n depends only on the vehicle tra- 
jectory. The exact values art' given in table If, 
although is given approximately by 

'>--[*■ (-"HT 

or 


0.S4 



The quantity 17, primarily a function of X, repre- 
sents the fractional reduction in mass loss due to 
latent heat. 

The relation between total mass loss and maxi- 
mum deceleration shows immediately the weight 
penalty incurred as the price of limiting the 
maximum deceleration to a low value, and it is 
concluded that the list* of a material with low 
ablation temperature is not appropriate* to vehicles 
which I lave high lift-drag ratios. It is seen in 
%ure 3 , for example, that the value of for 


L 

D 


0.5 is about twice* that 


for D- 


The effective heat capacity of the ablation 
material is writ ten 

(7?,. // ) x _ 0 (l — »,)-» Bi ii/ll, 

( f>. 2 

where (/A//) x=r , and t) are dim(*nsionless and do 
not vary appreciably with trajectory, as seen in 
figures 4 and 5 . Even when 17 ~() (negligible 
latent heat) tin* effective heat capacity of the 
material is 


l^OOa ■(//,„),_„ 

1 i>. 2 

when* 

0 . 3 5<^ ( II eff) x ^ 0 <T ♦ 5 


Again ( lioimporhmco of a!*- is soon ; when a—" —• , 

( (, p. 2 - 

for example, tin* effective heat capacity is be- 
tween 2,1)50 Bt ii/lh and 3,375 Btu/lb even when 
llu* latent heat is negligible. When this range of 
effective heat rapacities is compared with that 
for heat -sink metals of the order of 1,000 Btu/lb 
I he reduction in shield weight is quickly realized. 
Moreover, since the foregoing comparison does 
not depend on the ablation temperature, tin* 
advantage is enhanced when materials with low 
ablation temperatures are considered in view of 
the attendant reduction in insulation require- 
ments. 

For an ablation shield to perform successfully 
it must dispose of, rather than accumulate, 
heat energy. The prcahlution heating period 
should therefore be small compared with tin* total 
reentry time, and tin* heal accumulated should 
be a small fraction of that which would be ac- 
cumulated by a heat sink. Here, the deciding 
parameter is 


X-= 





1 1 has been shown that sublimation will not 
occur if x 2 >2 X 1 0 ft , that the ratio 


1 rcablation heating period 
Total heating period 

and that 

Heat accumulated by ablation shield 
Heat accumulated by heat-sink shield 


5(10 _fi )x* 

1 . 5 (l«- 3 )X 
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For the lifting vehicle sa > T ) Bie mass loss 

will be large if ablution is allowed to lake place 
during most of the reentry. Since sucli a vehicle 
would be cooled primarily by radiation, the use 
of a high-ablation-temperature material at the 
leading edge seems more appropriate. In such 
a design the ablation temperature should probably 
be near the imam radiation temperature of the 
vehicle, and the parameter x should be near the 
critical value 2X10® if ablation is to take place 
only near peak heating. The behavior ol the 
ablation material during this long preablation 
period may be of concern, however. 

When the heat conduction problem is consid- 
ered, an upper limit to the accumulation ol heat 
is found as 


Qt<- 2 v 2 J Btu/sq ft 

7T 

This^result is independent of the vehicle charac- 
teristics and heating experience except as they 
affect //. The temperature distribution satisfies 
the relation 



where the right-hand side is the limiting solution 
for negligible mass loss. From the foregoing 
expressions for Q and T it is seen that the amount 
of insulation will vary as // “. 

Thus, for a given vehicle the total shield weight 
required for sublimation and insulation varies 
approximately as the square root of the reentry 
time, or inversely as the square root of the maxi- 
mum deceleration. For ballistic vehicles ( 

5°) and manned capsules (0^ — (pi there- 
fore, the ablation shield offers an efficient way to 
dispose of heat continuously during reentry. 
For the lifting vehicle a high-ablation-tempera- 
ture material which would allow radiation from 
the surface for the greater part of reentry appears 
to be more appropriate; ablation would then 
take place for a limited time near the maximum 


heating condition or in case of an emergency 
maneuver. 


CONCLUDING REMARKS 


An approximate analysis has been made of abla- 
tion-shield requirements for reentry vehicles. I he 
type of shield considered was one of low ablation 
temperature and low thermal conductivity which 
produces no liquid film during ablation. It is 
shown that 

1. The total mass required for sublimation de- 
pends primarily on parameters which arc I unctions 
of the ratios 

Heating enthal py 

Gas shielding enthalpy 


and 


Shielding due to lat ent h eat- 
Gas-layer shielding 


2. For a given vehicle and shield the total 
mass loss varies as tin* square root ol the total 
time for reentry or inversely as the square root ol 
the maximum deceleration. 

;j. The heat accumulated is a small percentage 
of tha accumulated by a heat-sink shield, the 
percen age being determined by a single parame- 
ter winch combines the* effects ol the heating level 
experienced during reentry, the thermal capacity 
of the remaining shield, and the diflusivity ol the 
materi d. 

4. Tin 1 amount of insulation material also var- 
ies as he square root of the time or inversely as 
the square root of the maximum deceleration. 

From the foregoing dependence ol sublimation 
and ii solution requirements on deceleration and 
time of reentry it is concluded that the low- 
ablation-tcmpcrat uro shield should dispose of heat 
verv efficiently for lionlifting vehicles, but the 
limited use of a high-ablation-temperature shield 
at the leading edges is more appropriate for lilting 
vchicl *s, where the primary means ol cooling would 
he raciative. 


La nci. ;y Research Center, 

Nat (INAL Asia IN ATTICS ANI) Sl'ACE ADMINISTRATION, 
L no ley Field, Va., October 1 , thoiL 



APPENDIX 


FORMAL DEVELOPMENT OF THE DIFFERENTIAL EQUATIONS AND BOUNDARY CONDITIONS 

z = 0 is (ho sublimating surface. Thou the 


The boat conduction problem in a semi-infinite 
slab is represented by the following (‘([nation and 
boundary conditions. The heat conduction (‘([na- 
tion is 


Pt>(‘b 


57’ k 5*7’ 


(Al) 


dt " 5,: 2 

Tlu* initial conilition is 

7’=7’„ (/=<), //<<)) (A2) 

Tin* boundary conditions arc 

57’ 


and 


T-*T x 
T= T a 
d m 


a»/ 


(y=<>, t<f„) (Ada) 


( //- > 00 ; / <(/„) 


(» 


t>t 


qiO-h^ + L (H 

7 1 — >7’ 


Pr> 

m 


.) (A4a) 


(y=- m 'f>t u ) (A4b) 

(?/-->— 00 , tyt a ) (A4c) 

The heat -transfer rate </(/) is, in general, un- 
known and must be determined by an analysis of 
the aerodynamic heat -transfer problem. As stated 
in the body of the report the expression used 
here is 

'/(0 = '/ n (0 (*<0 (Aon) 


7(0 = 7„(0 ■ — <*c p ( T e — T a ) 


dm 

dt 


(tXa) (A5b) 


where r/ (j (/) is the aerodynamic heat-transfer rate 


to a nonablating body and ac if (T e — 


dm 


dt 

sents t ho shielding effect of the gas boundary layer. 


m 


The moving boundary ?/= — — is first elimi- 

Pb 


m 

nated by choosing a coordinate £=?/+ - so t ha t 

Pb 


(‘([nations ar(‘ made dimensionless by the intro- 
duction of new variables: 


T— T 


t 


' T„ r: T tr ? 

Equations (Al) to (A4) become 
ir \pJCb/ 


\Uf/ 


i/2 


(Ab) 


dr 

dr 


2 d tn dr d 2 r 
dt d£ = d£* 


= i) 
>0 


7o(t) J>r 

b( I a I co j V P Jr?, / d£ 

..... /<V/.Y /2 

'h(T a — 7 7 m ) ' k p tjcft,/ 


(A 7) 

(t = (), f^O) (ASa) 
(£ 00 , r>0) (ASb) 

(*=«>, (AO) 


5r 

5f 


rJ dm 
dt 


y' . Is-\-Qt(' p ( / I] i) / Cftfjr \'~ 

S Cl,{J a / oc ) \ Pbh'b ) 

0=». *>£) 


/■... 1 


0=°' '>?;) 


(A 10) 
(All) 


Fi 


tn 


’ dmidt is (‘(jiial to zero and the problem 

* f 

is the conventional one of finding the temperature 
distribution when the heat-transfer rate at the 

tn 


surface is given. 


For 


liowever, the differ- 


ential (‘([nation (A7) is nonlinear and there is an 
additional boundary condition corresponding to 

the additional unknown quantity ( ^ jf * 

The following approximate method is used to 
solve the foregoing system of (‘([nations. A mean 
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(%) % 


— e is used and the linear equa- 


tion which results from equation (A7) is solved 
to give r((.T,e), and e is then regarded as a function 
of r in I Ids solution. 

As a further simplifying approximation it is 
assumed that the preablation heating period is 
small compared with the total time of reentry; 

that is, -~<£l. Then y 1 at £ = 0 lor 0-^<A//, 

where / — 0 when 77— 0.995. 

. /Cft/rX 1 '" dm 

Equation (A7), with 6 replacing y - j J 


Ihm 


wit h 


Dr 
c>t 

r ■- 0 


dr dr d L V 

dT +( d! d? (Al2) 

1 ( 0 <r 1 ) (Alda) 

0 (| >- oo ,()<><!) ( A Idb) 


(t 0 , {<()) (Aide) 


The solution to equation (A12) is obtained by 
conventional Laplace transform methods and is 
written 


r <*' 


+J ( 2 T,, *“57'« ) ( Al4) 


The dimensionless heat-transfer rah 1 is obtained 
from equation (A 14) as 


( dr 

V of 


D,.. = 


1 

(pry ) 1 


d l 


e r * 
. (*rte • - 


(A 15) 


The third term e represents the dimensionless rate 
at which heat is absorbed by unit mass of material 
during its temperature list* from r - 0 to r- 1 and 
does not contribute 1 to the heat (*outent of the re- 
maining material. This heat content in dimen- 
sionless form is simply 


) ,T [(*t) -, '’ 2 < :- 4 t — ,* or ft* g r m 


dr 


0 


T ^ ^ 1 

( 1 “ T eric - r 1 


M- orf ; V'* 


When f— t-- 1, and €=e /? equation (A 14) 
reduces to 

>'f ( 1 orfc f/ 2 * )+' (Tfc (A17) 

which, a> has the limiting form 

»V=(Tfc(-|) (A1S) 

Tile asymptotic form for large e f is 

r r >e ( fl (A 19) 

Equatio l (Alb) for e ~ can be written 

Qf 2 /I 7T l/2 - Cr . 7T 1 « fA 

4 " J+2,; 2 ( ) 


“JTri? 

7T 


whore f >1 as and f- > - as e > co . 

^htf\ Ui /dm 


(A 20) 


The quantity e f - (^ is now replaced 


by 


/ c h tfY 2 n i 

f \phk,J t f 

i - u f 

W H - , /‘ (l-i») (-T) ( A21 > 

r a \pwb/ 


by using the relations 


( -V._Y" 

V/VlAv _ 


/a — 1 . 1 NA n _ m 


and 


777=0. K4 


' 4 (y) a-*) 


derived n tlie analysis of sublimation. 

Tin 1 product of dimensionless heat-transfer rate 
at the surface and the dimensionless heat content 
of the material (the quantity A which appeal's in 
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eq. (.Yli)) is found from i-qimtions (Al’>) and (Al(i) 
and is 


A — 



-m,. 



X 


[(If 


4 — - T (M'iV 


r 1/2 +- -vvtl 

e 




It mav bo verified that A--— when r- 0 tor all 

• 7T 

values of e and -1— >1 as t -■-><» tor all values of e. 
Thus it may bo expected that 

-^ 1^1 

7T 

That is, A is a slowly varying function of r. 
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